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Numerical simulations show that a long-range scalar interaction in a single species of massive dark 
matter particles causes voids between the concentrations of large galaxies to be more nearly empty, 
suppresses accretion of intergalactic matter onto galaxies at low redshift, and produces an early 
generation of dense dark matter halos. These three effects, in moderation, seem to be improvements 
over the ACDM model predictions for cosmic structure formation. Because the scalar interaction in 
this model has negligible effect on laboratory physics and the classical cosmological tests, it offers 
an observationally attractive example of cosmology with complicated physics in the dark sector, 
notably a large violation of the weak equivalence principle. 



PACS numbers: 98.80.-k, ll.25.-w, 95.35.+d, 98.65.Dx 

I. INTRODUCTION 

The ACDM model (cold dark matter with a cosmo- 
logical constant) passes demanding observational tests, 
but it may be possible to improve the agreement with 
observations by introducing an additional long-range in- 
teraction in the dark sector. In this paper, we consider 
a force law between two dark matter particles, each of 
mass to, arising from the potential 
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where (5 is of order unity and positive |36j| and the screen- 
ing length r s is on the order of 1 Mpc today and constant 
in co- moving coordinates. This type of force law, for par- 
ticles at a separation r much less than the Hubble scale, 
can arise from ordinary Newtonian gravity plus a mass- 
less scalar, where the scalar couples to the dark matter 
particles and to an additional relativistic particle species 
whose dynamics generates the screening length r s . The 
couplings of the scalar to the visible sector must be mini- 
mal in order not to conflict with laboratory tests of New- 
ton's law: the proposed modification Q applies only to 
dark matter particles. 

The idea of a long-range scalar interaction has been 
under discussion for a long time, in many contexts, as 
reviewed in 0, Q . The idea that the scalar interaction 
may couple only dark matter particles was introduced in 
Q . The field may relax to near minimum of its potential 
energy expressed as a sum over particle masses. In the 
example considered in Q this can eliminate the scalar 
interaction. In the scalar has a mass term with mini- 
mum offset from the minimum of the dark matter particle 
mass as a function of the scalar field. This produces a 
scalar interaction of the form (yi with r s a fixed proper 
length. The first analyses of the effect on dark mater halo 
formation and interactions are in 0, Q ■ In the model 
discussed here the scalar field finds its minimum poten- 



tial energy where the screening particles have close to 
zero mass, with the result that the cutoff length r s among 
the nonrelativistic dark matter species expands with the 
general expansion of the universe. This approach, which 
is developed in 0, H, , offers a variety of scenarios for 
cosmology: there may be several scalars with various cou- 
plings to several species of non-relativistic dark matter. 
For the purpose of a numerical exploration of what this 
approach can offer to physical cosmology, we focus here 
on a model with one non-relativistic dark matter species 
and one relativistic screening species. 

With /3 ~ 1, r s ~ 1 Mpc today, and r s constant not as a 
physical length but in co- moving coordinates, there is lit- 
tle effect on the classical cosmological tests, provided we 
assume that initial conditions are adiabatic and that the 
coupling to the screening particles suppresses evolution 
of the mean scalar field value. But our numerical simula- 
tions show that the effect on structure formation on the 
scale of galaxies is pronounced. The scalar force lowers 
the density of dark matter in the present-day voids be- 
tween the concentrations of large galaxies. It suppresses 
the rate of accretion of intergalactic debris onto galaxies 
at low redshift. And it increases the redshift of assembly 
of dark matter halos with comoving sizes smaller than r s . 

One could reproduce some of these arguably attractive 
results in a ACDM model with the initial mass distribu- 
tion contrived to match the present distribution of mass 
in our scalar interaction model. The match could be only 
approximate, of course, because the scalar force not only 
speeds up formation of galaxy-size halos but also also af- 
fects their interaction. And, in standard ACDM, it seems 
difficult to generate early haloes without also having sub- 
stantial late time accretion. 

In fjnjwe review the dark matter dynamics in the scalar 
field model under consideration and the motivation for 
the model from string theory. In ffllll we present the 
results of numerical simulations of structure formation. 
Our simulations have modest spatial resolution and num- 
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bers of particles, compared to what can be done, and 
should therefore be regarded as only a first exploration 
of what the model has to offer. Nevertheless, as discussed 
in section fllll CI we can conclude that the effects we do 
establish are observationally attractive. In ijlVI we com- 
ment on issues related to massive halos which we cannot 
analyze using our simulations, but which we hope will be 
explored in future work. 



II. DYNAMICS 

The form (JJJ is identical [^3] to the standard 
parametrization of hypothetical fifth force corrections to 
Newton's law of gravity in the visible sector (see for ex- 
ample the review article [Tpj') . The interesting regime for 
the strength parameter (3 is also the same as in discus- 
sions of a fifth force: (3 ~ 0(1). But the interesting scale 
for the screening length in our case is r s ~ 1 Mpc, vastly 
larger than the regime r s ~ 100 fim of interest in modern 
fifth force experiments. And, crucially, we assume that 
the scalar force acts only in the dark sector: visible mat- 
ter interacts with the dark matter only through ordinary 
gravity. 

It is essential that we assume that the self-interaction 
potential V{(f>) for the scalar field can be neglected: a 
significant non-zero V'(0) would drive cf> away from at 
late times, and non-zero V"(0) (with V'(0) = 0) amounts 
to a mass for <p which dominates over the screening ef- 
fect of the relativistic particles at late times. Non-zero 
higher derivatives of V at <f> = might endanger our 
story in subtler ways which have not been fully probed. 
These requirements on V(<f>) are in conflict with standard 
notions of field-theoretic naturalness. We nevertheless 
claim some motivation for our model from string theory 
and supersymmetry, as we briefly review in flll Bl 

Our further discussion in sections III Al and III Bl of the 
physics behind the force law Q includes some recapit- 
ulation of earlier discussions 0, H, El , presented here in 
the interest of a self-contained presentation. 



A. The screening mechanism 

Let one dark matter species have mass m — ycf>, while 
a second species has mass y s (f>. The couplings y and y s 
are positive dimcnsionless constants. Assume that the 
particles of mass y s (f> have a much larger number density, 
n <C n s , so that as the field moves to minimize the energy 
in particle masses it is pulled to < <f> <C m. That 
makes the screening particles relativistic and gives the 
dark matter particles masses that are close to m. The 
particle dynamics, independent of spin or statistics, can 
for present purposes be modeled as gasses of classical 
point-like particles with the action 
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where the sum is over all the particles, and m a (<f) = 
ui — y<f> or y s cj> depending on the species of particle a. 

During structure formation the scalar field dynamics 
can be treated in a quasi-stationary approximation: 



V 2 4> = 4>/r 2 s -yn{v,t) 
where V 2 is the spatial laplacian and 



r s = Vts/y 2 



(3) 
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The last term in © accounts for the non-relativistic 
particles in a hydrodynamic approximation. The pre- 
vious term, 4>/r 2 sl follows by noting |2| that the source 
term for <f> for a particle with speed v includes a factor 
ds/dt = \/l — v 2 , and that for quasi-static configurations 
of (j) the screening particle energy, e s = ra s j\J\ — v 2 = 
y a (f)/\/\ — v 2 , is nearly independent of position. Elimi- 
nation of yl — v 2 in favor of e s results in the screening 
length J3J|. The energy e s does change with time, scal- 
ing with the expansion of the universe as e s cx a(t) . 
The screening length thus scales as r s oc a(t), that is, the 
comoving length is constant. 

The scalar field produced by a single dark matter par- 
ticle at distance r <C r s is = y/Airr. The force this field 
exerts on another dark matter particle is the negative of 
the gradient of the mass m — y<p, that is, F = yV0, 
so we see that the particles are attracted with force 
F = y 2 /4irr 2 at r <C r s . This means the ratio of the 
scalar and gravitational forces of attraction of two dark 
matter particles is 



y 
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(5) 



The relations Q and JSJ summarize how the parameters 
of the potential Q emerge from the dynamics we have 
added to the dark sector. 

The effect of the scalar interaction on the evolution of 
the mass density contrast 5 = Sp/p in linear perturba- 
tion theory is simply expressed in terms of the Fourier 
amplitudes Sk(t). In the approximation that all the mass 
is in the dark matter, the evolution equation is 



6 k + 2-4 = AnGp 
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where k and r s are constant in comoving coordinates. In 
the Einstein-de Sitter model, modes with wavenumber k 
grow as 



«J k cx t p , 
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to be compared to the usual power law, p — 2/3, at kr s <ti 
1. The scalar interaction thus causes earlier development 
of small-scale structure. We comment on the possible 
implications in 3IVI 

The value of the interaction cutoff length r s is limited 
by the allowed energy density p s — e s n s in the relativistic 
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screening matter. The ratio to the mean mass density is 

Bl = hn(Hr s ) 2 S 2 , (8) 
P 2 

where the ratio of mean number densities of screening 
and dark matter particles is 

S= V -^. (9) 
yn 

To avoid violating the standard model for the origin 
of light elements at high rcdshift we must choose pa- 
rameters so the mass density in screening matter is less 
than that of about two low mass neutrino families [Tl| . 
or p s 0.5aT 4 , where the present background radiation 
temperature is T — 2.725 K. For the Hubble and matter- 
density parameters H = 70kms _1 Mpc _1 and f2 = 0.3, 
this condition in eq. iJSJ says 

(3S 2 r 2 s < 3000, (10) 

where r s is measured in megaparsecs. 

We can write the present screening particle energy as 

^- = ly^Smr,) 2 * 5 x lQ-V 2 Sr 2 s , (11) 

where M P1 = 1/V8ttG w 2.4 x 10 18 GeV is the reduced 
Planck mass. 

The growing concentrations of dark matter pull <f> away 
from zero. Where y s <p exceeds e s the screening particles 
are excluded — their mass would exceed their energy — 
and the scalar force is unscreened. The condition for this 
to happen is most simply stated in the limit where the 
characteristic width of the concentration is greater than 
r s , so V 2 is small compared to <f>/r 2 . Then eq. J3J may 
be rewritten as 

1 + 5 = S{^1 -v 2 ) < S, (12) 

where the matter density contrast is S = n/n — 1. When 
the length scale of the dark matter concentration is 
smaller than r s the screening matter is excluded from 
regions with a larger value of 6. Our numerical simula- 
tions ignore this exclusion of screening particles, so we 
are underestimating the the scalar attraction in regions 
where the dark matter mass density is large. Since the 
screening lengths we shall consider are greater than the 
sizes of the regions where the dark matter density is large 
enough to exclude the screening particles our simulations 
can ignore this somewhat complicated situation and fo- 
cus on a spatially uniform r s which is constant over time 
in comoving coordinates. 

B. String theory motivations 

In the string theory literature, there is an extended 
development, starting with [121 Il3| and continuing with 



substantial works on "brane-gas cosmology" [3^ , of pos- 
sible cosmological implications of strings or branes wrap- 
ping cycles in extra dimensions whose sizes and shapes 
are described by massless scalars in four dimensions. For 
the most part, this literature is concerned with effects 
that are uniform over the three spatial dimensions that 
we observe. For example, finite uniform number densities 
of strings with momentum or winding number around a 
circular fifth dimension tend to stabilize the size of the 
circle ll|, because the masses of the winding and mo- 
mentum strings are convex functions of the canonically 
normalized scalar controlling the size of the circle. 

The global picture emerging from this literature is 
the following. At tree level, most known compactifica- 
tions of string theory have a number of moduli (massless 
scalars) and a number of stable, heavy, point-like objects 
whose ten-dimensional description is in terms of wrapped 
or stretched branes (or, sometimes, Kaluza-Klein exci- 
tations). The masses of these objects are functions of 
the scalars, and these functions typically vanish each at 
a different point in moduli space. Sometimes (as for 
strings wrapping a circle) the vanishing point is at an 
infinite distance in moduli space. Tree-level exchanges 
of the massless scalars mediate long-range forces in four 
dimensions that are comparable to the force of gravity: 
in our terminology, j3 ~ 0(1). The reason is that these 
scalars are themselves gravitons (or perhaps superpart- 
ners of gravitons), but with polarization tensors oriented 
in the extra dimensions. There are also gauge inter- 
actions between the wrapped branes whose strength is 
roughly gravitational and whose ten-dimensional origins 
include gravitons whose polarization tensor has one in- 
dex in the compact directions and one in the four that 
we observe. But gauge interactions will be of little in- 
terest to us because their effect on structure formation 
is less pronounced than scalar interactions of comparable 
strength j^. 

This picture clearly presents a rich and fascinating 
field of possibilities for dark sector physics, particularly 
when we contemplate string compactifications with sev- 
eral dozen complex scalar moduli and dozens if not hun- 
dreds of stable brane configurations. But there are some 
potential problems too: 

1. After supersymmetry is broken, a potential is typ- 
ically generated for all the moduli on at least the 
meV scale. 

2. Even before supersymmetry breaking, there typ- 
ically are couplings of the scalars to the visible 
sector which would violate fifth-force constraints — 
unless said scalars have a mass greater than roughly 
an meV. 

3. If we want the dark matter to be wrapped branes, 
some mechanism must be specified to generate the 
right number density of them at early times. 

To implement a screening mechanism in a cosmologically 
interesting range of parameters (in particular, (3 ~ 1 and 
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r s ~ 1 Mpc) , the typical energy e s of a screening par- 
ticle should be on the order 10~ 7 y B M P \ today. Arrang- 
ing for such a large per-particle energy could be added to 
point 3 in the above list as a potential problem: ordinary 
production mechanisms of light particles seem unlikely 
to lead to large e s . 

In H, 0, some suggestions were made about how 
points 1 and 2 might be addressed. Before supersym- 
metry breaking, point 2 is mostly an issue of controlling 
the Kahler potential K{(f)\(j}) for the scalars which re- 
main massless. It is sufficient for the Kahler potential 
to have a minimum with respect to these scalars at the 
point in moduli space where we sit today. Then scalar 
couplings to the visible sector will start at dimension six, 
and estimates || show that solar system tests of New- 
ton's law of gravity are not compromised (and table-top 
laboratory tests would be entirely unaffected). It seems 
that we are asking rather a lot when we demand that 
the Kahler potential should be quadratic in the mod- 
uli just when one of the dark matter species becomes 
massless. But, as discussed in in certain supersym- 
metric theories with non-abelian gauge fields and ad- 
joint matter, the gauge invariance forbids linear terms 
in the Kahler potential along the Coulomb branch of 
the moduli space: precisely what we want. In such a 
situation, the massless scalar field </> (arising from the 
adjoint matter) would be complex, the screening parti- 
cles would include the non-abelian gauge bosons, and the 
whole ensemble could be arranged to be weakly coupled 
in the infrared by including enough matter fields (for in- 
stance, fundamentals and anti- fundamentals). It seems 
very plausible that some suitably constructed set of co- 
incident or intersecting branes could realize the physics 
outlined here, though an explicit example that includes a 
non-relativistic species with suitable couplings would be 
desirable. 

Clearly, the sticking point is supersymmetry breaking. 
Even if one assumes that the fundamental scale of super- 
symmetry breaking is Asusy = lOTeV and that its ef- 
fects are felt in the dark sector only through gravitational 
mediation, it is still difficult in conventionally understood 
scenarios to push scalar masses significantly below the 
gravitino mass, which is roughly 



^3/2 
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Without supersymmetry, the expectation from natural- 
ness is that scalar masses would be much bigger. 

The estimate (II 31) is the cleanest naturalness argument 
against our proposal. There are, however, some reasons 
to believe it can be circumvented. In the duality between 
strings in anti-de Sitter space (AdS) and a conformal field 
theory (CFT) on its boundary [lj, Il5l llq , the existence 
of Lorentz- invariant operators with dimension A ~ 0(1) 
in the CFT guarantees that there is a scalar field in AdS 
with mass 



v/A(A-4)/L. 



(14) 



where L is the radius of curvature of AdS — the analog 
of the Hubble scale. Supersymmetry is not invoked in 
any way in this argument, though the conclusion would 
come as a surprise to proponents of naturalness when L 
is much bigger than the string scale. It is conjectured 
|l7j | that similar arguments play out in the setting of an 
expanding universe entering a de Sitter phase, with the 
conclusion that there should be scalars with Hubble scale 
masses. The constraints on possible interactions of such 
scalars with dark matter are not well understood by the 
present authors, but they are presumed to be less severe 
than for Goldstone bosons. 

Also, there have been suggestions that scalar masses 
may be made extremely small if strong coupling renders 
mass terms irrelevant [l||i though there is again some 
question of what interactions are possible for such scalars. 

The need for highly energetic screening particles (point 
3 above) may prove to be the thorniest issue. Large e s 
might be arranged if the screening particles are the decay 
product of some very heavy and moderately long-lived 
particle — say a Planck mass particle that decays at z = 
lQ 7 /y s . It should also be recalled from [H, that one 
may assume more moderate values of e s and the mass m 
of the non-relativistic particles — bringing the problem of 
abundances perhaps more within reach of standard ideas 
like thermal production and freeze-out — if one lowers y s 
and y proportionally, so as to preserve f3 and r s . 

The simple model (|2J with one non-relativistic species 
and one relativistic one might arise straightforwardly as 
one characteristic behavior of many different underlying 
string constructions: the relativistic species corresponds 
to the wrapped branes that happen to be most numerous 
in the initial conditions; the non-relativistic species cor- 
responds to the next most numerous type of wrapped 
branes; and other massive species comprise negligible 
number and energy densities. 

In summary, if certain favorable features of classical 
string theory or unbroken supersymmetry are preserved 
at low energies in the dark sector, then the scalar in- 
teractions we study may arise from a broad variety of 
string compactification scenarios. Thus, it is just possi- 
ble that some simple ideas from string theory are closer 
to astrophysical observation than has previously been ap- 
preciated. 



III. NUMERICAL SIMULATIONS OF 
STRUCTURE FORMATION 



We turn now to exploring the consequences of the 
scalar interaction on the formation of structure in the 
expanding universe. This section deals with numerical 
simulations of structure formation on scales ~ 100 kpc 
to ~ 5 Mpc. We comment on the demanding issue of 
structure formation on smaller scales in WV\ 
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A. Numerical Methods 

The treatment of the particle motions given their pe- 
culiar accelerations is standard. To simplify notation in 
the computation of accelerations we write the scalar field 
as 

<f> = -AnGm<i)/y. (15) 

Then the force on a dark matter particle is 

F = -mV$ s + = -mV($ g + 0$), (16) 

where $ ff is the Newtonian gravitational potential and 
(3 is the ratio of the scalar to gravitational forces on 
scales small compared to r s (eq. |5])- In this notation the 
Fourier transforms of the gravitational and scalar field 
equations are 

fc 2 $ s (k) = (k 2 + r s " 2 )$(k) = -4irGpS(k), (17) 

where p = mn is the dark mass density, 8 — pj 1 p — 1 is 
the local density contrast. We have neglected the mass in 
baryons and assumed that the spatial average of V 2 cj> is 
zero. An FFT particle accelerator is easily adjusted from 
the standard computation of the gravitational potential 
Q g to compute the scalar field $. 

To get some indication of the baryon distribution we 
follow the motion of trace particles under the gravita- 
tional acceleration g = — V<I> 5 alone. These particles 
have the same number and initial conditions as the dark 
matter. We call these particles "baryons" though their 
motions do not take account of hydrodynamical forces. 

The particle motions were simulated with a particle- 
mesh (PM) FFT N-body code kindly provided by E. 
Bertschinger [l^. All simulations use 128 3 particles in 
a cubic box of 128 3 grid points for FFT computations. 
The initial conditions are represented by displacements 
from a cubic lattice to represent Gaussian density fluc- 
tuations with the ACDM power spectrum normalized so 
that the linear value of the rms density fluctuations in 
spheres of radius 8ft -1 Mpc is 0.9 at the present time. 
(The Hubble parameter is H = 100ft km s" 1 Mpc" 1 .) All 
runs started at redshift z = 100. The N-body code is 
run in practice with 17 = 0.3 and A = 0, but the output 
times are scaled to fi = 0.3 and A = 0.7 using the recipe 
in HJ. 

B. Simulations 

We present results from simulations with two values of 
the box widths, L — 10 and 50ft. -1 Mpc. For a given set 
of the initial conditions, simulations were run for several 
values of (3 and r s . For the larger box size L — SOh' 1 Mpc 
we obtained six simulations with different initial phases 
for each pair of (3 and r s . For L = 10ft" 1 Mpc we ran two 
simulations with different initial phases. 

The main theme in the following discussion is the effect 
of the scalar interaction in suppressing the amount of 
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FIG. 1: Present particle distributions in a slice 0.5/i _1 Mpc 
deep through a simulation box of width 50/i -1 Mpc. The ini- 
tial conditions are the same in each simulation. The panel 
at the top right shows the distribution in the absence of the 
scalar force. In the other simulations (3 = 1. The baryons in 
the lower right panel respond to the mass distribution in the 
lower left panel. 



dark matter debris between and around the massive dark 
matter halos. This is illustrated first in two maps of 
particle distributions and then in some statistics. The 
observational appeal is discussed in £1111 CI 

Figure ^ shows the effect of a scalar force of attrac- 
tion with /3 = 1 for r s = 1.56ft," 1 Mpc and 3.12ft" 1 Mpc 
in simulations with box size L = 50ft _1 Mpc. The con- 
trol case with no scalar force is shown in the slice in the 
upper right, and the dark matter distributions for the 
same slice with the same initial conditions and the two 
values of r s are shown in the left panels. One sees that 
the scalar interaction produces more prominent massive 
halos and leaves less dark matter between the halos and 
filaments. Since the particles representing baryons re- 
spond only to the gravitational force, it is not surprising 
that the baryon distribution in the lower right panel is 
less tightly clustered than the mass distribution in the 
lower left that produced it and more strongly clustered 
than the mass distribution in the control plot directly 
above it. 

Another aspect of the suppression of debris between 
massive halos is illustrated in figure|21 which shows exam- 
ples of the late time evolution of halos with and without 
the scalar interaction. For the purpose of our analysis, 
a "massive halo" is identified by the standard friends- 
of-friends (FOF) algorithm. In this algorithm, a par- 
ticle is in a group if it is within a prescribed linking 
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M=log[l+<5] 



FIG. 2: Illustrations of the suppression of accretion at low 
redshift. Each panel at z — shows the projected positions in 
a halo selected by the friends-of-friends linking algorithm with 
a linking length of 0.2 times the mean particle separation. The 
panels labeled z = 1 show the positions at that epoch of the 
particles selected in the linking list at z = 0. The comoving 
width of each panel is 3/i -1 Mpc and its height is 2/i -1 Mpc. 



FIG. 3: The distributions of the density contrasts in dark 
matter and baryons smoothed with a top-hat spherical win- 
dow of radius 1.5A -1 Mpc at the present epoch. The standard 
model is the solid curve, the dotted curve shows the effect 
of the scalar force with r s = 0.78/i _1 Mpc, and the dashed 
curve shows r s — 1.56/i -1 Mpc. The simulation box width is 
50ft' 1 Mpc. 



length of some other particle in the group. The sim- 
ulation box width for the identification of these halos 
is L = 50/i -1 Mpc, the mean interparticle separation is 
w = L/128, and the linking length in the FOF algorithm 
is 0.2w = 80/i _1 kpc. The comoving width of each panel 
is 3/i _1 Mpc and its height is 2/i _1 Mpc. The particle 
mass is 10 9 ' 7 /i _1 M©, and the halo masses in figure |21 are 
in the range 10 117 to 10 12 - 7 /i -1 M©. The panels labeled 
z = show the present positions of the particles iden- 
tified by the linking list, and the panels labeled z = 1 
show the positions of the same sets of particles at the 
earlier epoch. The examples without the scalar interac- 
tion show the prominent low redshift accretion predicted 
by the ACDM cosmology. Turning on scalar interactions 
leads to less accretion at low redshift, not because the 
scalar forces discourage accretion (the opposite is true) 
but because there is less debris to accrete. Loosely speak- 
ing, accretion finishes earlier when scalar interactions are 
turned on. 

The behavior illustrated in figures[T]and|21is quantified 
by the probability distribution function in the density 
contrast found within a randomly placed sphere, which 
we estimate from the density obtained by convolving the 
spatial particle distribution through a spherical window. 
The sphere radii are r = 1.5/i — 1 Mpc and r — hh^ 1 Mpc 
in the distributions in figures |3| and These sample 
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FIG. 4: The same as figure|3|for window radius 5/i -1 Mpc. 
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length scales are intermediate between the nominal halo 
sizes of galaxies, r ~ 300/i _1 kpc, and the typical 
sizes of voids, r ~ 15/i _1 Mpc. The mass distributions 
for the baryons, shown in the top panels, are not greatly 
affected by the scalar force, as one sees in figure ^ At 
both sphere radii the scalar force increases the probability 
of finding a sphere with large density contrast, but the 
much larger effect is the increased probability of finding 
a nearly empty sphere that contains density less than 3% 
of the mean. 

For a measure of structure formation on smaller scales 
we use a halo mass function n(M) defined such that 
n(M)d log 10 M is the mean number density of halos in the 
relevant mass range. Figure|3]shows n(M) for halos iden- 
tified in the L = lOh^ 1 Mpc simulations using the FOF 
algorithm. The linking parameter is 15/i _1 kpc and the 
particle mass is 10 7 ' 6 /i _1 .Mq. All halos with more than 
20 particles are included in the calculation of n(M). Fig- 
ure shows histograms of the mass function in the simu- 
lations with (3=1. The dotted and solid lines correspond 
to n(M) at z = and z = 1, respectively. The abundance 
of halos in the mass range 10 115 — f0 12 5 Af© is consis- 
tent with the observed abundance of galaxies thought to 
reside in halos in this mass range [21j. The halo mass 
function can be used to quantify the recent history of 
merging in halos with M ~ 1O 12 M : the scalar interac- 
tion shifts merging activities to higher masses relative to 
the standard dynamics without scalar interaction. This 
is encouraging as it mitigates some of the problems as- 
sociated with the recent intense merging predicted by 
standard ACDM halos of mass ~ 10 12 M Q (see fllTTty . 

The simulations with L = 10/i^ 1 Mpc are too small 
to assess the abundance of halos in the mass range cor- 
responding to rich galaxy clusters. Therefore we have 
also computed n(M) from the simulations with L = 
50/i Mpc. All of these larger simulations yield simi- 
lar prediction for the abundance of cluster halos. This is 
expected since even though r s is of the order of the virial 
radii of rich clusters, most of the matter in these objects 
collapsed only recently from much larger comoving dis- 
tances and therefore its dynamics has not been affected 
by the scalar interaction. 

In light of the limited resolution of our simulations, it 
is prudent to consider an alternative measure of the halos 
identified by our FOF algorithm. For the no scalar sim- 
ulation, figure IHI shows the numbers of particles N(<Rq) 
within radii Ro = 100/i _1 kpc and Ro = 300/i _1 kpc of 
the center of mass of each FOF halo. This figure gives 
evidence that the more massive FOF halos have the mass 
structures and density contrasts of conventional relaxed 
halos. The evidence is in two parts. First, at iVfof > 1000, 
we see that A^(<300/i" 1 kpc)/A(<I00/i" 1 kpc) » 3, in- 
dicating agreement with the commonly discussed form 
p ~ I/r 2 . Second, at iVf f = 1000 (corresponding to 
M = 1O 1Q -^ _1 M ) and R = 100/i _1 kpc, the density 
contrast is 5 ~ 100, which is close to the nominal density 
contrast at virial equilibrium. For Nf a i < 100, the signif- 
icant vertical scatter comes from small FOF halos which 
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FIG. 5: The halo mass function computed in a simulation 
box with width 10/i -1 Mpc. The dotted line is the present 
function, the solid line the mass function at redshift 2 = 1. 
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FIG. 6: Comparison of the particle number in a halo identified 
as a linking list with the number of particles within distances 
100 and 300/i." 1 kpc for the simulation with no scalar interac- 
tion. 



are close to big ones. In sum, the behavior illustrated in 
figure El is an encouraging indicator of the robustness of 
the notion of halos as identified by the FOF algorithm. 

Figure \7\ shows the cross- and auto-correlation func- 
tions of the present distributions of baryons and dark 
matter. The bend in the functions at separations 
r ^ 1 Mpc is an artifact of the force resolution in the sim- 
ulation, and should be disregarded, but we nevertheless 
expect that the effect of the scalar interaction indicated 
by the differences of the functions at 0.3 ,$ r 1 Mpc is 
meaningful. The most notable result is the absence of 
a feature at the length scale r s of the scalar force, de- 
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FIG. 7: The dark matter and "baryon" correlation functions. 



spite the scalar force of attraction at smaller scales. It is 
interesting also that at r ~ 300ftXkpc this measure of 
the baryon clustering is not much affected by the scalar 
interaction. And it is worth noting that the scalar in- 
teraction increases the slope of the mass autocorrelation 
function at 1 to 3ft. -1 Mpc. A tantalizing prospect is that 
this might reduce scale-dependent bias. 



C. Observational Situation 

At low redshifts, the scalar interaction lowers the den- 
sity of dark matter in voids and suppresses the rate of 
accretion of dark matter by galaxies. In this section we 
comment on how these two related effects may improve 
upon ACDM in matching observations. 

Voids between the concentrations of large galaxies con- 
tain plasma clouds with atomic hydrogen surface densi- 
ties ~ 10 13 cm -2 detected as Lya resonance absorption 
lines , but voids contain strikingly few isolated dwarf 
or irregular galaxies p3|. That is not what we might 
have expected from simulations of the dark mass distri- 
bution in the ACDM cosmology, as illustrated in figure 1 
of |24|. The analysis in [24| indicates that the distribu- 
tions of distances to the nearest low mass halo from giant 
and low mass halos are not very different, consistent with 
the observations. This is difficult to interpret, however, 
because the distances in the simulation are an order of 
magnitude larger than what is found for this statistic 
applied to galaxy catalogs [2^, and because our visual 
impression is that the simulation presented in pij shows 



distinctly more low mass halos between the concentra- 
tions of giants than are observed. This could be because 
the low mass halos in voids contain too little gas or stars 
to be observable. But that does not agree with the ob- 
servations of nearby dwarfs at ambient densities close to 
the cosmic mean. One would wonder why dwarfs in the 
voids, where the ACDM cosmology predicts densities just 
an order of magnitude lower than the mean, are not sim- 
ilarly visible. 

Perhaps the voids appear empty because a scalar force 
has pushed most of the dark matter out of them. Since 
the voids would have grown out of smaller density minima 
(as measured in co-moving coordinates), their growth 
would have been assisted by the scalar force even at the 
relatively short ranges r s we have considered, and indeed 
figure 01 shows an appreciable effect on the abundance of 
very low density regions with diameter lO/i" 1 Mpc. In 
this model it would not be surprising to see HI clouds 
in the voids, as one sees in figure ^ It wm require more 
detailed analyses to check whether the parameters can 
be chosen for consistency with, among other things, the 
constraints on the mass fraction in the baryons left in the 
voids and the observations of empty voids considerably 
larger than 10 Mpc. 

In the ACDM cosmology one expects considerable ac- 
cretion at low redshift. This is illustrated in figure 2 of 
[25). Figure |21 of the present paper shows that the scalar 
interaction distinctly suppresses accretion of dark mat- 
ter at low redshift. The observations seem to favor sup- 
pression; we mention three aspects. First, the evidence 
reviewed in j2(J, from the ages, chemical abundances, 
and spatial distributions of the stars in our Milky Way 
galaxy, is that this system has not been substantially 
disturbed b y a ccretion since redshift z ~ 1. Second, 
Blitz et al. [23 show that if the high velocity HI clouds 
were falling into the Milky Way from distances ~ 1 Mpc 
it would be a good match to the accretion expected in 
the standard ACDM cosmology. However, as discussed 
in |2jj, one does not observe HI emission around other 
galaxies — hence a challenge to ACDM if we accept the 
interpretation of [27]]. Third, in the judgement of Abadi 
et al. 29] the effect of accretion at low redshift on the 
thin discs of other spiral galaxies is "not grossly incon- 
sistent with current data" but "is worryingly difficult to 
accommodate within this general scenario." We conclude 
that although the standard cosmology is not ruled out 
by the observations, it does require that the ordinary- 
looking Milky Way galaxy is quite unusual, and that the 
recent accretion by other galaxies is well hidden. The 
more straightforward interpretation is that the ACDM 
cosmology might have to be adjusted, perhaps to include 
a scalar interaction in the dark sector. 

We have not explored whether the suppression of de- 
bris around massive halos is accompanied by a reduction 
in the numbers of dwarf satellites of a giant galaxy, which 
arguably would be observationally desirable [3(j ■ Check- 
ing this requires more detailed simulations. 
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IV. MASSIVE HALOS 

With one massive dark matter species the scalar inter- 
action must cause primeval mass fluctuations on scales 
;$ r s to grow into nonlinear massive halos distinctly 
sooner than in the standard cosmology. Our simulations 
are not adequate to show whether this is a challenge for 
the scalar model, or perhaps an advantage; we can only 
present some qualitative considerations and suggestions 
for further work. 

The growth of the first generation of halos commences 
when the mass density in matter becomes dominant, at 
redshift ~ z cq = 3500 for = 0.3 and h — 0.7. The mass 
density contrast subsequently grows as S ~ a(t) on scales 
larger than r s and, if f3 = 1, as S ~ a(t) 3 ^ 2 on smaller 
scales (eq. Q). This means the growth factor to z = 20 
in the standard model, when significant structure forma- 
tion commences, is accomplished in this scalar model at 
redshift z ~ 100. 

We get some understanding of the nature of this early 
generation of massive halos by noting that the Fourier 
amplitudes of the mass distribution at wavelengths larger 
than r s are not much affected by the scalar interaction. 
This is because the short range of the scalar force gathers 
dark matter from distances ~ r s , evacuating a compen- 
sating hole around a growing halo, so that the gravita- 
tional attraction of matter at larger distances is not much 
different from the standard model. The baryons would 
tend to lag behind the early growth of the dark mat- 
ter density contrast in these halos, because the baryons 
are attracted only by gravity. Thus we expect the for- 
mation of a class of gravitationally isolated, compact, 
baryon-poor halos. The mass autocorrelation functions 
in figure[7|have adequate resolution to be sensitive to the 
mass in the early halos. We expect that the absence of a 
feature at r s is a result of the growth of the clustering of 
matter gathered from larger scales at lower redshifts. 

More detailed simulations are needed to determine the 
dark matter mass fraction in these early halos and the ac- 
creted baryon mass under the assumption that nongrav- 
itational forces may be neglected. It will be exceedingly 
difficult, but important, to estimate the early star forma- 
tion rate in these halos, in order to decide whether it can 
be consistent with the reionization history indicated by 
the Wilkinson Microwave Anisotropy Probe [3l| . Also to 
be considered is whether the early generation of massive 
halos helps promote the strikingly early development of 



the black holes that power the SDSS z ~ 6 quasars |32l |. 



V. CONCLUDING REMARKS 

The numerical simulations presented here of a simple 
version of the scalar dark matter interactions considered 
in 0, [s. 9] open up a window to the non-linear regime. 
It is mainly in this regime, so crucial to comparison with 
observations, that the scalar force is expected to make a 
contribution to structure formation. 

We have shown that a dynamically screened scalar in- 
teraction offers clear relief to two quite troublesome as- 
pects of the ACDM cosmology. First, it promotes lower 
density of dark matter in the voids. This agrees with the 
observation that void dwarf and irregular galaxies are 
rare Second, it suppresses accretion of intergalactic 
debris onto galaxies at low redshifts. This agrees with 
the evidence that present-day galaxies by and large act 
like isolated island universes (as discussed, for example, 
in [23, |33| . And it is a fascinating possibility that the 
physics that improves these aspects of the otherwise very 
successful ACDM cosmology could be a window to string 
theory. 

Issues remain. If the scalar interaction preserves stan- 
dard estimates of the mass density as a function of ra- 
dius in the relaxed parts of a dark matter halo, then 
the earlier structure formation on scales r s would in- 
crease the halo mass function without affecting the issue 
of dark matter central cusps in galaxies (as discussed in 
pjof). But does it? The critical checks of the halo den- 
sity run and halo mass function require more extensive 
numerical studies. 



Acknowledgments 

We have benefitted from discussions with Neal Dalai 
and David Hogg and from comments from the referee. 
We thank Antonaldo Diaferio for providing us with his 
FOF code. SSG was supported in part by the Depart- 
ment of Energy under Grant No. DE-FG02-91ER40671, 
and by the Sloan Foundation. AN thanks the Prince- 
ton Institute for Advanced Study for its hospitality and 
support. 



[1] C. H. Brans (1996), gr-qc/9705069. 

[2] G. R. Farrar and P. J. E. Peebles, Astrophys. J. 604, 1 

(2004), astro-ph/0307316. 
[3] T. Damour, G. W. Gibbons, and C. Gundlach, Phys. 

Rev. Lett. 64, 123 (1990). 
[4] T. Damour and A. M. Polyakov, Nucl. Phys. B423, 532 

(1994), hep-th/9401069. 
[5] B.-A. Gradwohl and J. A. Frieman, Phys. Rev. Lett. 67, 



2926 (1991). 

[6] B.-A. Gradwohl and J. A. Frieman, Astrophys. J. 398, 
407 (1992). 

[7] J. A. Casas, J. Garcia-Bellido, and M. Quiros, Class. 

Quant. Grav. 9, 1371 (1992), hep-ph/9204213. 
[8] S. S. Gubser and P. J. E. Peebles (2004), hep-th/0402225. 
[9] S. S. Gubser and P. J. E. Peebles (2004), hep-th/0407097. 
[10] E. G. Adelberger, B. R. Heckel, and A. E. Nelson, Ann. 



10 



Rev. Nucl. Part. Sci. 53, 77 (2003), hep-ph/0307284. 
[11] R. H. Cyburt, B. D. Fields, K. A. Olive, and E. Skillman 

(2004), astro-ph/0408033. 
[12] J. Kripfganz and H. Perlt, Class. Quant. Grav. 5, 453 

(1988). 

[13] R. H. Brandenberger and C. Vafa, Nucl. Phys. B316, 
391 (1989). 

[14] J. M. Maldacena, Adv. Theor. Math. Phys. 2, 231 (1998), 

hep-th/9711200. 
[15] S. S. Gubser, I. R. Klebanov, and A. M. Polyakov, Phys. 

Lett. B428, 105 (1998), hep-th/9802109. 
[16] E. Witten, Adv. Theor. Math. Phys. 2, 253 (1998), hep- 

th/9802150. 

[17] S. Gukov, A. Neitzke and C. Vafa, private communica- 
tion. 

[18] S. Thomas, private communication. 

[19] E. Bertschinger and J. M. Gelb, Comp. Ph. 5, 164 (1991). 
[20] A. Nusser and J. M. Colberg, MNRAS 294, 457 (1998). 
[21] K. Adelberger, C. Steidel, M. Pettini, A. Shapley, 

N. Reddy, and D. Erb (2004), astro-ph/0410165. 
[22] J. T. Stocke, J. M. Shull, and S. V. Penton (2004), astro- 

ph/0407352. 

[23] P. J. E. Peebles, Astrophys. J. 557, 495 (2001), astro- 
ph/0101127. 

[24] H. Mathis and S. D. M. White, Mon. Not. Roy. Astron. 

Soc. 337, 1193 (2002), astro-ph/0201193. 
[25] L. Gao, A. Loeb, P. J. E. Peebles, S. D. M. White, 

and A. Jenkins, Astrophys. J. 614, 17 (2004), astro- 

ph/0312499. 

[26] G. Gilmore, R. F. G. Wyse, and J. E. Norris, Astrophys. 

J. 574, L39 (2002), astro-ph/0207106. 
[27] L. Blitz, D. N. Spergel, P. J. Teuben, D. Hartmann, 

and W. B. Burton, Astrophys. J. 514, 818 (1999), astro- 

ph/9803251. 

[28] D. J. Pisano, D. G. Barnes, B. K. Gibson, L. Staveley- 
Smith, K. C. Freeman, and V. A. Kilborn, Astrophys. J. 



610, L17 (2004). 
[29] M. G. Abadi, J. F. Navarro, M. Steinmetz, and V. R. 

Eke, Astrophys. J. 597, 21 (2003), astro-ph/0212282. 
[30] S. Kazantzidis et al., Astrophys. J. 608, 663 (2004), 

astro-ph/0312194. 
[31] C. Bennett et al., Astrophys. J. Suppl. 148, 97 (2003), 

astro-ph/0302208. 
[32] X.-H. Fan et al., Astron. J. 128, 515 (2004), astro- 

ph/0405138. 

[33] W. J. G. de Blok, M. A. Zwaan, M. Dijkstra, F. H. Briggs, 
and K. C. Freeman (2001), astro-ph/0111238. 

[34] R. Easther, B. R. Greene, M. G. Jackson, and D. Kabat, 
JCAP 0401, 006 (2004), hep-th/0307233. 

[35] S. Watson and R. Brandenberger, JHEP 03, 045 (2004), 
hep-th/0312097. 

[36] A gauge interaction between identical particles would 
give negative (3. 

[37] In the fifth force literature, the strength and range pa- 
rameters are usually denoted a and A. 

[38] Because our work overlaps only slightly with the existing 
work on brane-gas cosmolo gy t o date, we feel justified 
in citing two recent papers |34l 135^1 from this literature 
whose list of references can serve as a guide to further 
exploration. 

[39] Gauge interactions cause like charges to repel, whereas 
scalar interactions cause them to attract. So, in the linear 
regime, an isocurvature mode describing separation of 
positive and negative gauge charges grows less quickly 
than the adiabatic mode. This is in contrast to the growth 
caused by scalar interactions, which is faster than in 
standard ACDM. Likewise, in the non-linear regime, the 
main dynamical effect of a gauge interaction would be to 
maintain approximate local charge neutrality, whereas a 
scalar force actually adds to the gravitational tendency 
to increase the density of overdense regions. 



